Impact of somatic mutations and transcriptomic alterations on cancer aneuploidy.

Aneuploidy has been recognized as one of hallmark of tumorigenesis since the early 20th century. Recent developments in structural variation analysis in the human genome have revealed the diversity of aneuploidy in cancer. However, the effects of gene mutation and expression in tumors on aneuploidy remain poorly understood. Here, we performed whole exome analysis of over 5,000 Japanese cancer cases and investigated the impact of somatic mutations and gene expression alterations on aneuploidy. First, we evaluated tumor content and genomic alterations that could influence aneuploidy. Next, we compared the aneuploidy frequency in 18 cancer types and observed that TP53 mutations were associated with the aneuploidy on specific chromosomes in colorectal and gastric cancers. Finally, we used expression analysis to isolate pathways involved in aneuploidy accumulation from tumors without TP53 mutations. Chromosomal instability and cell cycle aberration were associated with aneuploidy in TP53 wild-type tumors, and 26 commonly upregulated genes were identified in aneuploidy-high solid tumors without TP53 mutations. Among them, two cancer-related genes (CENPA and PBK) were involved in aneuploidy. Our integrated analysis revealed that both TP53 mutations and transcriptomic alterations independent of somatic mutations affect aneuploidy accumulation. Our findings will facilitate further understanding of diverse aneuploidies in the tumorigenesis.

Comparative analysis of tumor content estimation methods based on simu- lated tumor samples identified their impact on somatic variant detection in cancer whole genome sequencing.

Whole genome sequencing (WGS) in cancer genomics has become widespread with recent technological innovations, and the amount and types of information obtained from WGS are increasing rapidly. Appropriate interpretation of results is becoming increasingly important in clinical applications. This study aimed to evaluate the accuracy of tumor content estimation and its impact on somatic variant detection, using 100 simulated tumor samples covering 10-100% tumor content constructed from the sequencing data of cell line models. Extensive analysis revealed that the estimation results varied among computational analytical methods. Notably, there was a large discrepancy in low tumor content (≤ 30%). The reproducibility decreased in cases wherein chromosome-scale copy number changes were observed in normal cells. The minimum tumor content required to detect somatic alterations was estimated to be 10-30%. Identification of whole genome doubling was achieved with the lowest tumor content, followed by single nucleotide variation/insertion or deletion, structural variation, and copy number variation. Tumor content had a significantly higher impact on the false negatives than the false positives in variant calls. Results should be interpreted cautiously for samples wherein tumor content is a concern. These results can form the basis of developing important guidelines for evaluating cancer WGS.

Whole-genome and Epigenomic Landscapes of Malignant Gastrointestinal Stromal Tumors Harboring KIT Exon 11 557-558 Deletion Mutations.

Gastrointestinal stromal tumors (GIST) with KIT exon 11 deletions involving in codons 557-558 (KIT Δ557-558) exhibit higher proliferation rates and shorter disease-free survival times compared with GISTs with other KIT exon 11 mutations. We analyzed 30 GIST cases and observed genomic instability and global DNA hypomethylation only in high-risk malignant GISTs with KIT Δ557-558. Whole-genome sequencing revealed that the high-risk malignant GISTs with KIT Δ557-558 (12 cases) had more structural variations (SV), single-nucleotide variants, and insertions and deletions compared with the low-risk, less malignant GISTs with KIT Δ557-558 (six cases) and the high-risk (six cases) or low-risk (6 cases) GISTs with other KIT exon 11 mutations. The malignant GISTs with KIT Δ557-558 showed higher frequency and significance in copy number (CN) reduction on chromosome arms 9p and 22q, and 50% of them had LOH or CN-dependent expression reduction in CDKN2A. In addition, SVs with driver potential were detected in 75% of them, in which AKT3 and MGMT were recurrently identified. Genome-wide DNA methylation and gene expression analyses showed global intergenic DNA hypomethylation, SNAI2 upregulation, and higher expression signatures, including p53 inactivation and chromosomal instability, as characteristics of malignant GISTs with KIT Δ557-558 that distinguished them from other GISTs. These genomic and epigenomic profiling results revealed that KIT Δ557-558 mutations are associated with increased genomic instability in malignant GISTs. Significance: We present genomic and epigenomic insights into the malignant progression of GISTs with KIT exon 11 deletions involving in 557-558, demonstrating their unique chromosomal instability and global intergenic DNA hypomethylation.

Comparison of genetic susceptibility to lung adenocarcinoma and squamous cell carcinoma in Japanese patients using a novel panel for cancer-related drug-metabolizing enzyme genes.

The differences in genetic susceptibility to lung adenocarcinoma and squamous cell carcinoma remain unclear. We developed a customized, targeted gene sequencing panel for efficient and sensitive identification of germline variants, including whole-gene deletion types for cancer-related drug-metabolizing enzyme genes in lung adenocarcinoma and squamous cell carcinoma. The minor allele frequencies of the variants, confirmed as clinically significant in the Japanese population, did not differ significantly from those of normal participants listed in the public database. Genotype analysis comparing lung adenocarcinoma (n = 559) and squamous cell carcinoma (n = 151) indicated that the variants of DPYD (rs190771411, Fisher's exact test, P = 0.045; rs200562975, P = 0.045) and ALDH2 (rs568781254, P = 0.032) were associated with an increased risk of squamous cell carcinoma compared to adenocarcinoma. Conversely, whole-gene deletion of CYP2A6 was associated with adenocarcinoma but not squamous cell carcinoma. Notably, whole-gene deletion of CYP2A6 was confirmed in 22 patients with lung adenocarcinoma but not in any patients with squamous cell carcinoma. Most patients with whole-gene deletion of CYP2A6 were female non-smokers. The discovery of a whole-gene deletion of CYP2A6 in patients with lung adenocarcinoma may have an important role in clinical practice and advance our understanding of CYP2A6 germline variants and their association with carcinogenesis or their susceptibility to lung adenocarcinoma.

JCGA: the Japanese version of the Cancer Genome Atlas and its contribution to the interpretation of gene alterations detected in clinical cancer genome sequencing.

With the emergence of next-generation sequencing (NGS)-based cancer gene panel tests in routine oncological practice in Japan, an easily interpretable cancer genome database of Japanese patients in which mutational profiles are unaffected by racial differences is needed to improve the interpretation of the detected gene alterations. Considering this, we constructed the first Japanese cancer genome database, called the Japanese version of the Cancer Genome Atlas (JCGA), which includes multiple tumor types. The database includes whole-exome sequencing data from 4907 surgically resected primary tumor samples obtained from 4753 Japanese patients with cancer and graphically provides genome information on 460 cancer-associated genes, including the 336 genes that are included in two NGS-based cancer gene panel tests approved by the Pharmaceuticals and Medical Devices Agency. Moreover, most of the contents of this database are written in Japanese; this not only helps physicians explain the results of NGS-based cancer gene panel tests but also enables patients and their families to obtain further information regarding the detected gene alterations.

BMP4 and PHLDA1 are plausible drug-targetable candidate genes for KRAS G12A-, G12D-, and G12V-driven colorectal cancer.

Despite the frequent detection of KRAS driver mutations in patients with colorectal cancer (CRC), no effective treatments that target mutant KRAS proteins have been introduced into clinical practice. In this study, we identified potential effector molecules, based on differences in gene expression between CRC patients carrying wild-type KRAS (n = 390) and those carrying KRAS mutations in codon 12 (n = 240). CRC patients with wild-type KRAS harboring mutations in HRAS, NRAS, PIK3CA, PIK3CD, PIK3CG, RALGDS, BRAF, or ARAF were excluded from the analysis. At least 11 promising candidate molecules showed greater than two-fold change between the KRAS G12 mutant and wild-type and had a Benjamini-Hochberg-adjusted P value of less than 1E-08, evidence of significantly differential expression between these two groups. Among these 11 genes examined in cell lines transfected with KRAS G12 mutants, BMP4, PHLDA1, and GJB5 showed significantly higher expression level in KRAS G12A, G12D, and G12V transfected cells than in the wild-type transfected cells. We expect that this study will lead to the development of novel treatments that target signaling molecules functioning with KRAS G12-driven CRC.

Genomic profiling of multiple tissues in two patients with multiple endocrine neoplasia type 1.

Multiple endocrine neoplasia type 1 (MEN1) is an autosomal dominant tumor syndrome. This hereditary cancer is caused by germline variants in MEN1. Two patients with MEN1 were identified via whole exome sequencing and gene expression profile analysis, conducted for 5,063 patients with various types of cancers. We obtained multiple tumors from each patient; tumors derived from these two MEN1 patients had a loss of the normal MEN1 allele and frequently chromosomal copy number changes. Thus, we investigated whether structural variants were present in the MEN1 patient genomes. Whole-genome sequencing revealed no catastrophic rearrangements, and the tumor samples had very low somatic variants. The two patients had germline variants in MEN1 and some chromosomal copy number changes including on chromosome 11. The only pathogenic variant detected was the MEN1 germline variant, and chromosomal rearrangements led to tumorigenesis in somatic cells. Furthermore, the MEN1 tumor samples displayed a specific signature characterized by T:A>C:G transition. Studies of multiple tumors obtained from single patients are rare in hereditary cancer syndromes, and our results provide insights that the second hit of the tumor suppressor gene MEN1 may be caused by a gross genome rearrangement, not a small insertion and deletion, nor a change in epigenetic regulation.

Mutational concordance analysis provides supportive information for double cancer diagnosis.

BACKGROUND: Mutation analysis using next-generation sequencing highlights the features of tumors with somatic alterations. However, the mutation profile of double cancer remains unclear. Here, we analyzed tumors derived from the same patient using whole exome sequencing (WES) to investigate the coherence of somatic mutations in double cancer. METHODS: First, the tumor mutational burden (TMB) was investigated using WES of 5521 tumor specimens from a Japanese pan-cancer cohort. The frequencies of mutation concordance were then compared in these cancers. Finally, we calculated the expected value of mutational concordance fitting a Poisson distribution to determine the relationship between double and metastatic cancers. RESULTS: In all, 44, 58, and 121 paired samples were diagnosed as double cancer, multifocal lesions (derived from identical tissues), and metastasis, respectively. Our analysis revealed that common somatic mutations were almost entirely absent in double cancer, whereas primary tumors and metastatic foci harbored several identical alterations. Concordance of the mutation profile in the same patient reflects the tumor origin and development, suggesting the potential for identifying double cancer based on common somatic mutations. Furthermore, according to a Poisson distribution, double cancer could be discriminated based on paired samples from the same patient. The probability of double cancer with more than 10 mutations was ≤1 part-per-billion (ppb, 10- 9). In multifocal lesions, 74% of tumor pairs accumulated ≤10 common mutations, implying a difference in tumor origin within identical tissues. CONCLUSIONS: These findings indicate that counting common somatic mutations can indicate the differences in origin between tumors derived from the same patient. Our mutation coherence analysis can thus provide beneficial information for diagnosing double cancer.

Characterization of tumors with ultralow tumor mutational burden in Japanese cancer patients.

Tumor mutational burden analysis using whole-exome sequencing highlights features of tumors with various mutations or known driver alterations. Cancers with few changes in the exon regions have unclear characteristics, even though low-mutated tumors are often detected in pan-cancer analysis. In the present study, we analyzed tumors with low tumor mutational burden listed in the Japanese version of The Cancer Genome Atlas, a data set of 5020 primary solid tumors. Our analysis revealed that detection rates of known driver mutations and copy number variation were decreased in samples with tumor mutational burden below 1.0 (ultralow tumor), compared with those in samples with low tumor mutational burden (≤5 mutations/Mb). This trend was also observed in The Cancer Genome Atlas data set. In the ultralow tumor mutational burden tumors, expression analysis showed decreased TP53 inactivation and chromosomal instability. TP53 inactivation frequently correlated with PI3K/mTOR-related gene expression, implying suppression of the PI3K/mTOR pathway in ultralow tumor mutational burden tumors. In common with mutational burden, the T cell-inflamed gene expression profiling signature was a potential marker for prediction of an immune checkpoint inhibitor response, and some ultralow tumor mutational burden tumor populations highly expressed this signature. Our analysis focused on how these tumors could provide insight into tumors with low somatic alteration that are difficult to detect solely using whole-exome sequencing.

Japanese version of The Cancer Genome Atlas, JCGA, established using fresh frozen tumors obtained from 5143 cancer patients.

This study aimed to establish the Japanese Cancer Genome Atlas (JCGA) using data from fresh frozen tumor tissues obtained from 5143 Japanese cancer patients, including those with colorectal cancer (31.6%), lung cancer (16.5%), gastric cancer (10.8%) and other cancers (41.1%). The results are part of a single-center study called "High-tech Omics-based Patient Evaluation" or "Project HOPE" conducted at the Shizuoka Cancer Center, Japan. All DNA samples and most RNA samples were analyzed using whole-exome sequencing, cancer gene panel sequencing, fusion gene panel sequencing and microarray gene expression profiling, and the results were annotated using an analysis pipeline termed "Shizuoka Multi-omics Analysis Protocol" developed in-house. Somatic driver alterations were identified in 72.2% of samples in 362 genes (average, 2.3 driver events per sample). Actionable information on drugs that is applicable in the current clinical setting was associated with 11.3% of samples. When including those drugs that are used for investigative purposes, actionable information was assigned to 55.0% of samples. Germline analysis revealed pathogenic mutations in hereditary cancer genes in 9.2% of samples, among which 12.2% were confirmed as pathogenic mutations by confirmatory test. Pathogenic mutations associated with non-cancerous hereditary diseases were detected in 0.4% of samples. Tumor mutation burden (TMB) analysis revealed 5.4% of samples as having the hypermutator phenotype (TMB ≥ 20). Clonal hematopoiesis was observed in 8.4% of samples. Thus, the JCGA dataset and the analytical procedures constitute a fundamental resource for genomic medicine for Japanese cancer patients.

Driver gene alterations and activated signaling pathways toward malignant progression of gastrointestinal stromal tumors.

Mutually exclusive KIT and PDGFRA mutations are considered to be the earliest events in gastrointestinal stromal tumors (GIST), but insufficient for their malignant progression. Herein, we aimed to identify driver genes and signaling pathways relevant to GIST progression. We investigated genetic profiles of 707 driver genes, including mutations, gene fusions, copy number gain or loss, and gene expression for 65 clinical specimens of surgically dissected GIST, consisting of six metastatic tumors and 59 primary tumors from stomach, small intestine, rectum, and esophagus. Genetic alterations included oncogenic mutations and amplification-dependent expression enhancement for oncogenes (OG), and loss of heterozygosity (LOH) and expression reduction for tumor suppressor genes (TSG). We assigned activated OG and inactivated TSG to 27 signaling pathways, the activation of which was compared between malignant GIST (metastasis and high-risk GIST) and less malignant GIST (low- and very low-risk GIST). Integrative molecular profiling indicated that a greater incidence of genetic alterations of driver genes was detected in malignant GIST (96%, 22 of 23) than in less malignant GIST (73%, 24 of 33). Malignant GIST samples groups showed mutations, LOH, and aberrant expression dominantly in driver genes associated with signaling pathways of PI3K (PIK3CA, AKT1, and PTEN) and the cell cycle (RB1, CDK4, and CDKN1B). Additionally, we identified potential PI3K-related genes, the expression of which was upregulated (SNAI1 and TPX2) or downregulated (BANK1) in malignant GIST. Based on our observations, we propose that inhibition of PI3K pathway signals might potentially be an effective therapeutic strategy against malignant progression of GIST.

Mutational burden and signatures in 4000 Japanese cancers provide insights into tumorigenesis and response to therapy.

Tumor mutational burden (TMB) and mutational signatures reflect the process of mutation accumulation in cancer. However, the significance of these emerging characteristics remains unclear. In the present study, we used whole-exome sequencing to analyze the TMB and mutational signature in solid tumors of 4046 Japanese patients. Eight predominant signatures-microsatellite instability, smoking, POLE, APOBEC, UV, mismatch repair, double-strand break repair, and Signature 16-were observed in tumors with TMB higher than 1.0 mutation/Mb, whereas POLE and UV signatures only showed moderate correlation with TMB, suggesting the extensive accumulation of mutations due to defective POLE and UV exposure. The contribution ratio of Signature 16, which is associated with hepatocellular carcinoma in drinkers, was increased in hypopharynx cancer. Tumors with predominant microsatellite instability signature were potential candidates for treatment with immune checkpoint inhibitors such as pembrolizumab and were found in 2.8% of cases. Furthermore, based on microarray analysis, tumors with predominant signatures were classified into 2 subgroups depending on the expression of immune-related genes reflecting differences in the immune context of the tumor microenvironment. Tumor subpopulations differing in the content of infiltrating immune cells might respond differently to immunotherapeutics. An understanding of cancer characteristics based on TMB and mutational signatures could provide new insights into mutation-driven tumorigenesis.

Tumor mutational burden analysis of 2,000 Japanese cancer genomes using whole exome and targeted gene panel sequencing.

Tumor mutational burden (TMB) is an emerging characteristic in cancer and has been associated with microsatellite instability, defective DNA replication/repair, and response to PD-1 and PD-L1 blockade immunotherapy. When estimating TMB, targeted panel sequencing is performed using a few hundred genes; however, a comparison of TMB results obtained with this platform and with whole exome sequencing (WES) has not been performed for various cancer types. In the present study, we compared TMB results using the above two platforms in 2,908 solid tumors that were obtained from Japanese patients. For next-generation sequencing, we used fresh-frozen tissue specimens. The Ion Proton System was employed to detect somatic mutations in the coding genome and to sequence an available cancer panel that targeted 409 genes. We then selected 2,040 samples with sufficient tumor cellularity for TMB analysis. In tumors with TMB-high (TMB ≥ 20 mutations/Mb), TMB derived from WES correlated well with the estimated TMB (eTMB) based on panel sequencing, whereas TMB in the remaining tumors showed a weak correlation. In particular, eTMB was overestimated in tumors with low-frequency mutations, resulting in the accumulation of EGFR mutations not being discriminated as a feature of lung cancer with low-frequency mutations. The eTMB in tumors harboring POLE mutations and microsatellite instability was not overestimated, suggesting that panel sequencing could accurately estimate TMB in tumors with high-frequency mutations such as hypermutator tumors. These results may provide helpful information for interpreting TMB results based on clinical sequencing using a targeted gene panel.

A novel MLH1 intronic variant in a young Japanese patient with Lynch syndrome.

Lynch syndrome, an autosomal dominantly inherited disease, is characterized by an increased risk of developing colorectal cancer. We found a novel germline variant of MLH1 (IVS6+2T>C) that caused Lynch syndrome in a young Japanese patient who had multiple colorectal cancers. Accurate diagnosis will be highly beneficial in clinical practice for surveillance and genetic counseling of patients and their relatives.

Germline and somatic genetic changes in multicentric tumors obtained from a patient with multiple endocrine neoplasia type 1.

Multiple endocrine neoplasia type 1 (MEN1) is a hereditary cancer syndrome caused by germline mutations of the MEN1 gene located in chromosome 11q13. In patients with MEN1, multicentric tumors develop in the involved organs; however, precise evaluation of genetic changes in these multicentric tumors has not been performed. In the present study, using whole-exome sequencing, we analyzed germline and somatic genetic changes in blood cells, two pancreatic endocrine tumors and one duodenal tumor obtained from a patient with MEN1 gastrinoma. We found that this patient possessed a novel germline mutation of the MEN1 gene [NM_137099.2:c.1505dupA (p.Lys502Lysfs); the localization was Chr11:64572134 on Assembly GRCh37], in which an adenine insertion in codon 502 of the MEN1 gene resulted in a frame shift and a premature stop codon. In terms of heterozygosity, the mutated allele was heterozygous in blood cells, hemizygous in the two pancreatic tumors and homozygous in the duodenal tumor. Immunohistochemical staining confirmed that only truncated menin protein accumulated in the nucleus of the tumor tissues. Further evaluation of tumor-specific somatic mutations in two pancreatic tumors did not detect single-nucleotide variations (SNVs) in 609 cancer-associated genes designated by the COSMIC cancer gene census, suggesting that the germline MEN1 mutation and resultant loss of heterozygosity played a major role in tumorigenesis. In the duodenal tumor, in addition to the germline MEN1 mutation, single-nucleotide variations in two cancer-associated genes were found. Further studies are required to clarify the role of these somatic single-nucleotide variations in the progression of MEN1 tumors.

Comprehensive characterization of genes associated with the TP53 signal transduction pathway in various tumors.

The TP53 signal transduction pathway is an attractive target for cancer treatments. In this study, we conducted a comprehensive molecular evaluation of 907 patients with cancer in Japan to identify genomic alterations in the TP53 pathway. TP53 mutations were frequently detected in many cancers, except melanoma, thymic tumors, gastrointestinal stromal tumors, and renal cancers. The frequencies of non-synonymous single nucleotide variants (SNVs) in the TP53 family members TP63 and TP73 were relatively low, although genes with increased frequencies of SNVs were as follows: PTEN (11.7%) in breast cancer, CDKN2A (11.1 and 9.6%) in pancreas and head and neck cancers, and ATM (18.0 and 11.1%) in liver and esophageal cancers. MDM2 expression was decreased or increased in patients with mutant or wild-type TP53, respectively. CDKN1A expression was increased with mutant TP53 in head and neck cancers. Moreover, TP63 overexpression was characteristically observed in squamous cell carcinomas of the lung, esophagus, and head and neck region. Additionally, overexpression of TP63 and TP73 was frequently observed in thymomas. Our results reveal a spectrum of genomic alterations in the TP53 pathway that is characteristic of many tumor types, and these data may be useful in the trials of targeted therapies.

Whole exome sequencing detects variants of genes that mediate response to anticancer drugs.

Certain interindividual differences affecting the efficacy of drug treatment and adverse drug reactions are caused by genetic variants, and their phenotypic effects differ among ethnic groups. In this study, we used whole exome sequencing (WES) systematically to identify germline mutations that influence the activities of drug-metabolizing enzymes, as well as that of a transporter. We analyzed DNA isolated from blood samples from 2,042 Japanese patients with diverse cancers. We identified sequence variants of CYP2B6 (rs3745274), CYP2C9 (rs1057910), CYP2C19 (rs4986893), CYP2C19 (rs4244285), TPMT (rs1142345), NAT2 (rs1799930), NAT2 (rs1799931), UGT1A1 (rs4148323), COMT (rs4680), ABCB1 (rs1045642), and CDA (rs60369023). Wider application of WES will help to determine the effects of mutations on the activities of proteins encoded by drug response genes, and the information gained will accelerate the development of personalized therapies for patients with cancer. Moreover, this knowledge may provide clues for preventing cancer before the onset of symptoms.

Optimizing an ion semiconductor sequencing data analysis method to identify somatic mutations in the genomes of cancer cells in clinical tissue samples.

Identification of causal genomic alterations is an indispensable step in the implementation of personalized cancer medicine. Analytical methods play a central role in identifying such changes because of the vast amount of data produced by next generation sequencer. Most analytical techniques are designed for the Illumina platform and are therefore suboptimal for analyzing datasets generated by whole exome sequencing (WES) using the Ion Proton System. Accurate identification of somatic mutations requires the characterization of platform-dependent error profiles and genomic properties that affect the accuracy of sequence data as well as platform-oriented optimization of the pipeline. Therefore, we used the Ion Proton System to perform WES of DNAs isolated from tumor and matched control tissues of 1,058 patients with cancer who were treated at the Shizuoka Cancer Center Hospital. Among the initially identified candidate somatic single-nucleotide variants (SNVs), 10,279 were validated by manual inspection of the WES data followed by Sanger sequencing. These validated SNVs were used as an objective standard to determine an optimum cutoff value to improve the pipeline. Using this optimized pipeline analysis, 189,381 SNVs were identified in 1,101 samples. The analytical technique presented here is a useful resource for conducting clinical WES, particularly using semiconductor-based sequencing technology.

Integrated next-generation sequencing analysis of whole exome and 409 cancer-related genes.

The use of next-generation sequencing (NGS) techniques to analyze the genomes of cancer cells has identified numerous genomic alterations, including single-base substitutions, small insertions and deletions, amplification, recombination, and epigenetic modifications. NGS contributes to the clinical management of patients as well as new discoveries that identify the mechanisms of tumorigenesis. Moreover, analysis of gene panels targeting actionable mutations enhances efforts to optimize the selection of chemotherapeutic regimens. However, whole genome sequencing takes several days and costs at least $10,000, depending on sequence coverage. Therefore, laboratories with relatively limited resources must employ a more economical approach. For this purpose, we conducted an integrated nucleotide sequence analysis of a panel of 409-cancer related genes (409-CRG) combined with whole exome sequencing (WES). Analysis of the 409-CRG panel detected low-frequency variants with high sensitivity, and WES identified moderate and high frequency somatic variants as well as germline variants.

Next generation sequencing approach for detecting 491 fusion genes from human cancer.

Next-generation DNA sequencing (NGS) of the genomes of cancer cells is contributing to new discoveries that illuminate the mechanisms of tumorigenesis. To this end, the International Cancer Genome Consortium and The Cancer Genome Atlas are investigating novel alterations of genes that will define the pathways and mechanisms of the development and growth of cancers. These efforts contribute to the development of innovative pharmaceuticals as well as to the introduction of genome sequencing as a component of personalized medicine. In particular, chromosomal translocations that fuse coding sequences serve as important pharmaceutical targets and diagnostic markers given their association with tumorigenesis. Although increasing numbers of fusion genes are being discovered using NGS, the methodology used to identify such fusion genes is complicated, expensive, and requires relatively large samples. Here, to address these problems, we describe the design and development of a panel of 491 fusion genes that performed well in the analysis of cultured human cancer cell lines and 600 clinical tumor specimens.